import sympy as sp
from sympy import factorial, binomial, gamma, zeta, primepi, primerange

# 组合数学函数
print("10的阶乘:", factorial(10)) # 10的阶乘: 3628800
print("二项式系数 C(10, 3):", binomial(10, 3)) # 二项式系数 C(10, 3): 120

# 特殊数学函数
print("Gamma函数 Γ(5):", gamma(5)) # Gamma函数 Γ(5): 24
print("Riemann zeta函数 ζ(2):", zeta(2).evalf()) # Riemann zeta函数 ζ(2): 1.64493406684823

# 数论函数
print("小于100的质数个数:", primepi(100)) # 小于100的质数个数: 25
print("前10个质数:", list(primerange(1, 30))) # 前10个质数: [2, 3, 5, 7, 11, 13, 17, 19, 23, 29]
print("17是否是质数:", sp.isprime(17)) # 17是否是质数: True